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Abstract: The one loop effective potential in thermal gluodynamics has stable minima 
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1 Introduction 

The existence of a plasma phase of QCD has been confirmed by experiments at RHIC and 
LHC. Lattice simulations have been since long precious information about the equilibrium 
properties of the plasma. 

One question of interest is how the Stefan-Boltzmann gas of quarks and gluons at very 
high temperature develops into the strongly coupled plasma at lower temperatures, still 
above the critical temperature Tc- 

One can make the question more precise by looking at the free energy of the plasma 
as function of the order parameter, the Polyakov loop (Wilson line). This order parameter 
measures the surplus free energy of a heavy fundamental quark and reads for SU (2) 

^TrVex^p^J dtAo{x,t)^ = ^TtP{Ao{x)) (1.1) 

For a constant background Aq = 2TTqT'^ and the Wilson line is parametrized by 

^TrP(Ao(f)) = cos(7rg). (1.2) 
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For the sake of discussion we will exclude quarks from the system and do the calculation 
for gauge group SU{2). In that case the free energy is Z{2) symmetric, i.e. periodic mod 
1 and invariant under exchange of q and 1 — q. This symmetry is spontaneously broken at 
high temperature T, where the order parameter does not vanish. The one loop fluctuations 
around the translation invariant state Eq. (1.2) generate the free energy 

47r2 

ftr = fsB +-\q\f- (1-3) 

The minima of this free energy represent the actual equation of state of the plasma. 
At very high temperature this is the Stefan-Boltzmann (SB) law with the Polyakov loop 
being ±1. At temperatures just above the critical temperature the minima will become 
nearly degenerate (first order transition) or coalesce at g = | (second order transition, as 
indeed realized in SU{2), and where the Polyakov loop (1.2) vanishes.). 

Perturbative corrections to the potential do not change the location of the minima as 
is expected. One may ask: does the location change when semi-classical saddle points like 
those of thermal instantons or calorons are taken into account? Do these saddle points 
initiate the movement of the SB minima? 

To obtain the answer one has to compute the determinant of fluctuations around the 
saddle points in question. Its contribution to the free energy can be written on general 
grounds in terms of the running coupling g{T) as: 

ftot — ftr + fcal 

Ui = -2T^ (^^0 + nig + • • • ^ (1.4) 

The last factor in parenthesis is the linearized instanton density integrated over the 
caloron size. The dots are terms that render fcai Z{2) symmetric. 

To leading order in the coupling no,i are just numerical constants. They pick up a 
temperature dependence only in next to leading order, which warrants the calculation of 
two loop effects. 

The coefficient no was calculated in the seminal paper by Gross, Pisarski and Yaffe 
[1] (GPY in the sequel). It increases the Stefan Boltzmann pressure. The coefficient ni is 
computed in this paper. It has a negative sign, such that the free energy fcai increases. 
Therefore the minimum of the SB ground state at g = 0 stays stable. 

To further motivate the work done below we mention the work Diakonov and cowork¬ 
ers [2-5] (from now on DGPS). They have computed the caloron determinant as an ex¬ 
pansion in the inverse size ri 2 of the caloron. So their result is most reliable at large size 
ri 2 . Now the coefficients no,i are both obtained by an integral over the size. This integral 
converges in the ultra-violet due to asymptotic freedom (suppressing the small ri 2 contri¬ 
bution by a power), and in the infrared due to thermal screening. The latter suppresses 
exponentially the large size contribution, where the expansion of DGPS is most precise. 

In this paper we avoid the large size expansion and obtain an analytic result for the 
linearized density in Section 4, Eq. (4.22). 
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Figure 1. The potential (1.3) and in the blow up the correction due to (1.4) with a positive (dashed) 
and negative (dotted) linear term qni. The dotted line is realized by the caloron gas as shown in 
this paper. 

As shown below this coefficient can be expressed in terms of moments of the GPY 
result for the instanton density at g = 0. Recently the GPY result has been recast in 
analytic form [6], so ni is a very well controlled quantity. 

Yet another motivation for computing only the linear term niq comes from the depen¬ 
dence of the screening on the holonomy q. It is known since long [7, 8] that the screening 
mass in a constant background TrqTr^ calculated for diagonal fields is proportional to the 
second derivative of ftr 

(1.5) 

As screening is effective at large distance where the caloron by definition equals the holon¬ 
omy q, and the caloron fields are diagonal at large distance one would expect the screening 
to be of this type. This is indeed what DGPS do find for the screening of the caloron. 

This screening turns into anti-screening in a window centered around q = ^, where 
the potential is concave (see Fig. 1). Hence the contribution of the caloron will diverge in 
that window^. So one may limit oneself just as well to only the linear term, which by itself 
gives important information on the physics. 

Knowledge of that term determines the stability of the Stefan-Boltzmann minima. 
This is the strategy we will follow. 

Recapitulating, on the basis of our analytic result for this linear term we conclude that 
the minima do not change due to caloron fluctuations, they stay stable. Our statement is 
true to leading order in the running coupling g{T). 

feature common in mean field approximations. 


- 3 - 









The lay-out of the paper is as follows. In Section 2 the effective potential is introduced, 
and the main motivation for the paper explained. In particular the relevance of the term 
linear in a holonomy expansion is emphasized. In Section 3 the classical caloron solution 
is presented together with a discussion of the variation of the caloron with respect to the 
holonomy. Appendix A does furnish details. 

Section 4 contains the central part of our paper. In Section 5 our results are summarized 
and we compare them to those of DGPS, as far as possible. 

2 Effective potential and its behavior away from the Stefan- Boltzmann 
limit 

In this section we discuss the effective potential to one loop order in a constant Polyakov 
loop background cos( 7 rg). More precisely the effective potential is defined as a Euclidean 
path integral with the gauge invariant constraint that only configurations with 

dxTrP{Ao){x) = cos{'Kq) (2.1) 

are admitted [9] where obviously q is constant in space. The three volume V is supposed 
to be macroscopically large. An obvious choice for the vector potential is 

= 5^,o2vrTg^. (2.2) 

The classical action vanishes for all q, so is degenerate. The fluctuations around this 
translation invariant background do lift this degeneracy. 

The result to one loop [1, 10] and two loop [11-13] is a very simply ^-independent 
renormalization of the one loop result 

+ (2.3) 

The 0{g‘^) term is due to the two loop free energy diagrams in the background q 
The Stefan-Boltzmann free energy is 

fsB = ~T\ (2.4) 

Clearly the minima are in (7 = 0,1 where the loop takes the Z{2) values ±1, and this 
is still true to two loop order^. 

The question is now: can we find other local minima of the action with fluctuations 
furnishing for small holonomy a linear term in q? If we choose the caloron minimum then 
ftr changes into: 

ftr + fcal + • • • (2-5) 

^And a renormalized Polyakov loop insertion into the one loop free energy, which renders the result (1.4) 
gauge independent and kills all linear terms in q 

®This simple multiplicative renormalization is true for any group, with the factor 10 replaced by 5C2, 
the adjoint second order Casimir invariant[12] 
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with fcai as in Eq. (1.4). 

Then a negative coefficient ni will destabilize the minimum at (7 = 0. A positive sign 
stabilizes the minimum. See Fig. 1 and the blow up. 

Local minima of that kind are provided by (anti)-self-dual configurations with topo¬ 
logical charge ±1, and as boundary condition at spatial infinity that the Polyakov loop 
is cos( 7 rg). Such configurations have classical action , independent of the value of 

q. However one loop fluctuations lift the degeneracy, and indeed provide us with a linear 
term. However we find the sign to be positive and hence the minimum stays where it is: 
at (7 = 0 . 


3 The caloron with non-trivial holonomy and the Harrington-Shephard 
thermal instanton 


The caloron is an (anti)-self-dual solution of the Euclidean equations of motion, with the 
boundary condition that at spatial infinity the Wilson line P{Ao)(x) takes the value: 

^TrP(Ao)(f) = cos( 7 rg). (3.1) 

In the trivial holonomy case (g = 0), the caloron gauge field in the singular gauge is given 
by the Harrington-Shepard (HS) solution [14]: 


\ = -r]%d^^ogU 

H = 1 ri2f 
^ 1 sinh(r) 

r cosh(r) — cos(t) 


(3.2) 

(3.3) 

(3.4) 


where for convenience r and t have been rescaled r —)• 2'KrT and t —>■ 27rtT The’t Hooft 

symbol is given by: 

= V + (3-5) 

where = (1, —ir), f the Pauli matrices. 

In the general case ((7 7 ^ 0), the caloron can be seen as composite by two constituent 
monopoles [15-17]. At large separation ri 2 one is approximately a BPS [18, 19] monopole, 
the other is a twisted BPS monopole (in short a Kaluza-Klein (KK) monopole [20]). They 
have opposite magnetic charges and topological charges q and I — q respectively. So the 
solution has topological charge 1 and vanishing total magnetic charge. 

The action equals and is degenerate in the zero mode parameters and the 

holonomy q. 

Using rotational and translation symmetry we can put the two monopoles on the z axis 
and the centre of mass at the origin. The separation between the constituent monopoles 
is fixed by ri 2 : 


qXi -k (1 - q)X2 = 0 


X 2 — Xi = (0, 0, ri2) 

^The partial derivative will always be in the unsealed coordinates 


(3.6) 

(3.7) 
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Figure 2. Caloron with size ri 2 larger than the charge radii q ^ and g ^ of HS and KK monopoles 


The sizes of the BPS and KK monopoles are determined by the inverse of q and 1 — q 
respectively, so we have two different scenarios depending on whether <C ri 2 or ;§> 
ri 2 . In Figures 2 and 3 we present an example of each case: 

• Figure 2 corresponds to qr 12 ;§> 1: in this case both constituent monopoles are well 
separated. This is the case analysed in [2-5] with a large ri 2 expansion. 

• Figure 3 corresponds to qri 2 <C 1: in this case the KK monopole is embedded in the 
core of the BPS monopole. It is instructive to see how the rotational symmetry is still 
broken through terms of 0{q), i.e. as long as the core of the BPS monopole is finite, 
but large. Only when q = 0 the core is delocalized, and rotational invariance is exact. 
This physics is better reflected in a parametrization [21, 22] who use prepotentials 
that are specifically connected to the amount of breaking of rotational invariance. It 
reads 


{5abVs + eaSbVA + SsaSsbVs) - 


(3.8) 


This expression is quite useful since is a covariant derivative common to the scalar 
potentials Vs,A, 3 - Details are to found in Appendix A, where also the relation to the 
parametrization of KvB is given. Note that the tensor is S'0(3) invariant, the 
other two tensors only 50(2) invariant. Terms linear in q in the effective potential 
will be entirely controlled by the potential with the largest symmetry. Vs- The reason 
is that the trace in the loop calculation obliterates the terms with lower symmetry 
in V 5 and V 3 . An accidental symmetry cause V 3 to have vanishing 0{q) terms. 

The result in the limit qri 2 <C 1 for the caloron gauge potential can be written as 
(see Appendix A): 

(log (1 + ri 2 (l + ruq)/)) + g(NRI terms) + 0{q^) (3.9) 

The first term in Eq. (3.9), corresponds to a HS caloron with size ri 2 changed into 
^’ 12(1 + q'f' 12 )- The terms that are not 3d rotational invariant are indicate by NRI. 

This requires the computation of the variation of the caloron configuration with 
respect to q. Eq. (3.9) is the result of two observations: 




dqA‘l,{ri2,q)\q=o = -ruV] 


+ (NRI terms) 


(3.10) 

(3.11) 
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Figure 3. Caloron with size ri 2 smaller than the charge radius of the BPS monopole. 

The first relation is simply expressing the linearity of the HS caloron in terms of ri 2 , 
Eq. (3.4). The second relation is the result of the calculations in Appendix A. 

The two relations combine into 


5 qA“(ri2,g)|q=o = rl 2 dru^l{ri 2 ,q = 0) + (NRI terms), (3.12) 

which is nothing but the term linear in q in Eq. (3.9), using the HS potential Eq. (3.4). 
This is our key result. 

As we will see in Section 4 the NRI terms are washed out in the one loop integration. 
This means the linear term is determined, up to some minor modification (see Section 4) 
by the ri 2 derivative of the HS determinant, something we had already computed in our 
previous paper [6]. 

4 One loop contribution of the caloron 

The one loop contribution is given by the fluctuation determinant, integrated over the zero 
modes. Once the zero modes are factored out the fluctuation determinant of the spin one 
particle (written with a prime below) equals the square of that of the spin zero particle. 
The latter has no zero modes. Taking the ghost into account we get: 

Zi = j dCJiO exp det-^/2^-T);,^)det(-T)|;,) 

= J dCJ{0det-\-D^). (4.1) 
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The logarithmic divergencies are cured by normalizing Z\ with the determinant 
of the zero temperature instanton[23]^. 

This determinant in turn contains the one loop short distance effects for the coupling 

871-2 22 

' = — log(x) + log log terms. (4.2) 

g^{x) 3 

With A the Pauli-Villars scale equal to one gets. 

1 / 2 TtT 9 9 9 /X 

C stands for the eight zero modes parameters and J(C) is the zero mode measure 

J{C) = ^ + ri2q{l-q) (4.4) 

The covariant isospin one d’Alembertian for a scalar particle reads: 

{d^Sac + iA^eadc){dij,Scb + iA^Cceb) — ^ab^ ( 4 ^-^) 

where the potential is the caloron potential, Eqs. (A.l) and (A.2). The zero modes do 
not appear in the result of the determinant, except for the scale factor g. The zero mode 
integration then takes the form from Eq. (4.4): 

JdCJiC) = ^^J^ ^{'^ + n2qil-q))= J ^^{l + ri2q{l-q)) (4.6) 

ly is the volume of the SU{2) instanton. The space time volume V/T is due to the 
space time zero mode integration. 

The transition amplitude Zi is usually written as a product of the translation invariant 
part, and a part due to the caloron 

Zi = ZtrZi (4.7) 

= if) ifk)) “P (^9^) i ('‘■8) 

Here ^ g^^T) ) (~pfr)) q-iT / K) does represent the density per unit space 

time volume of calorons with size ri 2 , temperature T and holonomy q. To leading order in 
the coupling there is no T dependence in n{ri 2 , q,T/A), only in the prefactors. 

The one loop contribution of the translation invariant minima, i.e. with = 5^fi2TTqT'^ 
is ftr, as given in Eq. (2.3) 

Zi, = exp(-E/r/i,) (4.9) 

This bulk free energy gets a correction due to the caloron gas: 

^N+ 

^ N ' jv' =^^P(2^i) (4-10) 

N± 

®This normalization is used in Eq. (4.20) 
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Combining Eqs. (4.8), (4.9) and (4.10) one obtains for the total free energy 

ftot = ftr-2Zil{VlT). (4.11) 

So the free energy of the caloron becomes in terms of the density: 

{~9^) L 

l-OO 2 

n{ri 2 ,q,T/k) = C I ciri 2 r )^2 x (1 + ^2^(1 - ?)) exp(-logdet(-D^)) . (4.12) 

Jo 

The constant C equals 

C' = 2"ii/7r^. (4.13) 

The latter two equations tell us that the contribution to the bulk free energy equals 
minus twice the total density of calorons as a function of q, as defined in Eq. 4.8. Erom 
its definition in Eqs. 4.1 and 4.6 (positive caloron measure) and Eq. 4.8 one expects this 
total density to be positive, for all values of the holonomy q, for which it is defined. 

Eor vanishing holonomy one retrieves the density of Harrington Shepard calorons 
and the result of GPY[1]. The caloron gas produces a positive correction to the Stefan- 
Boltzmann pressure —fsB- 

4.1 The linear term in caloron effective potential and 3d rotational invariance 

Clearly the caloron free energy, Eq. 4.12, has a linear term in q. Eirst off all it receives 
a contribution from the measure term ~ ri 2 g(l — q), involving only the HS determinant. 
This contribution is clearly negative, so destabilizes the minimum. 

That the determinant exp(—log det(—11^)) delivers a linear term as well is shown by 
computing the derivative and setting q = 0: 

dgexp(-logdet(-T>^))|q=o = -dg (logdet(-T>^)) exp(-logdet(-i:>^)|q=o) 

= Tr [{dgD^){-D‘^)~^\g=o] exp(-logdet(-T»^)|g=o)(4.14) 

The trace stands for space integration, integration over the time period and tracing of 
color. 

The variation of the d’Alembertian at vanishing holonomy involves 

dgD^ = dqA^D^ + D^dqA^ (4.15) 

The last two equations show that the covariant derivative the propagator 

and the determinant are needed, but only at vanishing holonomy, i.e. they are rotational 

invariant. 

Note that the variation of the d’Alembertian on the right hand side of (4.14) is pro¬ 
jected on the rotational invariant HS propagator at g = 0®. In that projection only the 

®The HS propagator is only used as a regularization device, at coinciding points where it rotationally 
invariant. 
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rotational invariant part of dqA^, Eq. (3.12), survives. In other words, the linear term in 
Eq. (4.14) becomes, upon substitution of Eq. (3.12) and using Eq. (4.15), 

Tr = rf^Tv [{dr,,A,)D, +D^idr,,A,){-DY^)]^^^ + Tr[NRl] 

= r^^Tr 

The second equality follows because in the first line also the vector potentials A^ are zero 
holonomy potentials. The second set of terms vanishes for the mentioned symmetry reasons. 
The NRI terms are made explicit in Appendix A, Eq. (A.24). It is then straightforward 
to check they do vanish in the trace Eq. (4.14). 

So we conclude that the variation in q at zero holonomy is just the variation in ri 2 of 
the HS determinant, apart from an overall factor rf 2 - 

This means that the trace involved in Eq. (4.14), defining the coefficient of the linear 
term, is just as convergent as the derivative of the HS determinant with respect to the size 

n2^- 

4.2 The linearized instanton density 

To leading order the density n{ri 2 ,q,T/A) is temperature independent. The instanton 
density is positive for any value of A, and q, 

n(A, q, T) = no(A, g = 0, T) + (n(A, q, T) - n(A, O' = 0, T)) > 0. (4.17) 

In the preceding sections we have approximated the term in parenthesis by the linear 
term, 

n(A, q, T) = no(A, g = 0, T) + gni(A, T) H- (4.18) 

and noticed (see Eq. (4.16)) that we need only the HS determinant [6] 


det(—H 


- 2 ^ 


iq=o = exp ( - ( 16A + -ri 2 




R = 0.0129 1 + 0.899 


+ R 

(^)- 3 /. 


(4.19) 


R is uniformly quite small compared to the logarithm. 

The function A is the result of a logarithmically diverging A subtracted with the one 
instanton contribution Aq. Explicitly: 


A{ri2) 


f _ f 4 (0>.n.)^)b 

12164,2 ni T.nj ) 


with Ho the one instanton prepotential: 


Ho = 1 + 


7 


This convergence excludes terms like gloglg) in the effective action. 


(4.20) 


(4.21) 
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Some care has to be taken in applying (4.14) to the HS result. The reason is that 
Aq is not depending on the holonomy q. Hence it drops out when we vary the regulated 
determinant with respect to q, as in the first factor in (4.14). 

Hence using the input of (4.19) the integrated caloron density Zi becomes 


no+qni = C 


1^1 -k qri2{l - ri29ri2 -k ^ log(ri2) -k 

^ g-(l6A+fri2)_ 




(4.22) 

The term 


\ log(ri2) 

(4.23) 


results from omitting the one instanton subtraction. 


5 Results and comparison to other work 

Substituting Eq. (4.19) into Eq. (4.22) gives then the final result for the density ni(ri 2 ) 
and is shown in Eig. 4 together with the density no(ri 2 ). 

The outcome for the sign of ni is determined by a competition between the measure 
factor Eq. (4.6) and the contribution from Eq. (4.12) as is clear from Fig. (4) or from 
inspection of the integrand multiplying q in Eq. (4.22). It is the latter that wins, with the 
result 

dri2(no(ri2) + qni{ri2) = C(2.83854 - 17.5653g) (5.1) 

Substituting this in Eq. (1.4) we get for the total free energy: 

hot = hr - 2T^ exp Cl (^2.83854 - 17.5653g + • • • 

A 7]-^ 

hr = hr = fsB + - \q\?. (5.2) 

We see the caloron determinant gives a positive correction to the pressure —fsB and a 
positive linear term to the potential. Clearly the latter does not change when varying the 
temperature. 

Let us compare this to the work of DGPS. First in reference [4] there is a q log(g) term 
in caloron contribution. Such a term would give a logarithmically diverging coefficient 
ni- However we have seen that that coefficient is just as convergent as uq. On the other 
hand DGPS (see the relevant figure in [4]) do find from their numerical approach that 
below a certain temperature the caloron fluctuations do indeed destabilize the perturbative 
minima, but not above. That is, their first derivative of the effective potential at the origin 
is positive. So in that sense their result agrees with ours at very high temperature. 

Since the slope in Eq. (5.2) is positive for all values of T, the DGPS result that the 
slope at lower T switches sign must involve non-leading 0{g^) effects. 

The non-leading terms are partly due to two loop renormalization effects of the cou¬ 
plings. This is what is taken into account in the DGPS approach, are easy to compute and 
give indeed a negative 0{g^) term, and may indeed lead to a flip in the sign of Eq. (1.4) if 
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Figure 4. Upper panel: the HS caloron density no(ri 2 ) ; lower panel: the density of the coefficient 
ni{ri2) 

the coupling is large enough. However the two loop contribution (as yet to be calculated) 
may change that conclusion. 

The nice observation of DGPS that the screening term in logdet(—equals 

4ri2 Q - g(l - (5.3) 

is strictly valid when the two monopoles are well separated. To illustrate this remember 
our linear term gets its contribution from the region where the monopole with core 1/g is 
embedding the KK monopole, the one that survives in the HS limit. This is shown in Fig. 
3. In this regime we found, Eq. (3.9), that including the linear term just changes the size 
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of the HS caloron: 


ri2^ ri2{l + qri2)- (5.4) 

So starting from the HS screening |ri 2 we find a positive linear term That means 

switching on q increases the screening! This is precisely opposite to what the screening 
Eq. (5.3) tells us. This is another illustration of the difference between the two regimes, 
well separated monopoles and one monopole nested inside the other. DGPS devote an 
interesting effort to the connection between the two regimes. 

Finally we offer an amusing speculation. The potential increases from the SB minimum 
due to the caloron gas. However for larger q the potential may bend and go down creating 
a stable minimum at some temperature so that the SB minimum becomes meta¬ 

stable. The tunneling to this new minimum would give a glitch in the thermodynamic bulk 
functions. 

What might be the moral that switching on the caloron gas fails to initiate the insta- 
bilty of the Stefan-Boltzmann gas ? A cynic would say, that calorons are after all instantons 
and the latter are known since long to fail to give a first step to confinement. Proponents 
would say, that the monopoles inside the instanton might have given precisely that first 
step. Our computation shows that, alas, monopoles are not sufficient for destabilization. 
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A KvBLL caloron and its Taylor expansion around the HS caloron 

In this appendix we give some details about the determination of the Taylor expansion 
of the KvBLL caloron around the HS caloron in terms of the holonomy, q. In order to 
simplify this calculation, we found the notation introduced in [21, 22] very useful. Of 
course we compare with the result obtained from the Kraan-van Baal form of the solution. 
In Subsection A.l we remind the reader the KvB formula and the equivalence with the our 
notation. Finally in Subsection A.2 we explicit the computation of the Taylor expansion. 

A.l The Kraan - van Baal (KvB) solution explicit 

The solution in terms of KvB prepotentials becomes, with ri 2 = 2'k‘^T‘^q^^ qi = q,q 2 = 1 —9 
and rescaled coordinates [15, 17]: 
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■j sinh(giri) sinh(g 2 ?' 2 ) 1 / 2 rir 2 (cosh(( 7 iri) cosh(q 2 ?’ 2 ) — cost) 2 2 2 A 

rir 2 2\ sinh(giri) sinh(g 2 r 2 ) +^ 1+^2 ’’ 12 J 

sinh(giri) sinh(g2?’2) 


= ip + 

ip 


nr 2 


-ri 2 iri 2 + ri cosh(giri) + r 2 cosh(g 2 ?’ 2 )) 


(p = — (reduces to HS prepotential 11 in the —)• 0 limit.) 

-ip 

There is also a second -complex- prepotential: 


X = 


_ ri 2 sinh(giri) sinh(g 2 r 2 ) f iq^t _ n 


Ip 


rir2 


sinh(giri] 


+ e 


-iq 2 t_ 


r2 


sinh(g 2 ?' 2 ) 


The vector potential in terms of these prepotentials reads (r^ = 4(r^ ± 


(A.l) 


(A.2) 


= -vf,uT^du log (p - (p{T+ii^^d^x + r r]+^d^x*)- 


(A.3) 


Rotational invariance is gone and only rotations around the z-axis leave the caloron dipole 
configuration invariant. In particular the scalar potential equals 


Ao = vLdulogcp + (p{T^d-x + r d+x*) ■ 


(A.4) 


We wrote = 9i ± id 2 

The solution is invariant under q -n 1 — q and moreover t ■n —t for x- So the 
contribution to the free energy is Z(2) invariant. 

It is useful to rewrite the potential in terms of new prepotentials: 

= -n^FT^u {6abVs + ea3bVA + SsaSsbVs) - (A.5) 

When computing the variation in q it turns out that ^p starts out with a quadratic term in 
q. The other prepotentials have a complicated linear term, but the result in terms of Vs, 
Va, V 3 is relatively simple. It is given at the end of next subsection, Eq. (A.24). 


A.2 An insightful derivation of the variation 

In the notation introduced in [21, 22]. The SU(2) caloron gauge field is written just in 
terms of a 2 x 2 real matrix, V, which inverse is additive in constituents monopoles: 


\ = ^F{1, in) a^alduV 



+ h.c. 


(A.6) 


where 

y-i _ (Q.osh{qara) - cos{qat) (-l)“sin(gat) 

^^^4T^‘^smh.{mara) y (-1)“ sin(gat) cosh{qara) + cos{qat) 

and F-i = 1 - ^Tr(E). 
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Reordering Eq. A.6, the vector potential can be written as a covariant derivative, up 
to a common factor: 


{ 6 abVs + e^sbVA + iv 


where 


and 


V. = d.-d.logV, V = detE-\ F = 


ri2 F 
dvr^ V 


- (^■')22 = 


Tg COs{qat) 

27r2 sinh(gar’a) 
(-l)Va sin(gfat) 


Va-2(V-%,--^ 2»^smhfer.) 

U = 2(W')„ = g + J^ 


^■12 , ra(cOsh(ga?’a) “ COs{qat)) 


27r2 sinh(gara) 


where 


qi= q, 92 = 1 - 9 
rl = \x-Xa\^ 

This notation is related to the KvB formululation by: 

Ip = 

F 

= Att'^KV 

X = "yiVs + iVA) 


(A.7) 

(A.8) 

(A.9) 

(A.IO) 

(ATI) 


(A.12) 

(A.13) 


(A.14) 

(A.15) 

(A.16) 


where K = ^i°h( 9 in)sinh(g 2 r- 2 ) 

2 rir 2 

The main advantage of this notation is that F is already order q, so we only need the 
leading order of the other quantities. We start computing the Taylor expansion of F: 


F = 27r^9(n - 1) ( 1 - ri29 


r coth(r) 
ri2 


+ n^h ] ] + 0{q^) 


(A.17) 


where ns = ^ and h is defined in Eq. B.6. Since F is already order g, we need also the 
order due to the ^ term in Eq. A.9. In the case of V we only need the leading order: 


V = 


?'12 


1 n rcoth(r) 11 n^h , , 

+ 1 +--7 + 1 + 0{q) 


dvr^ Vri29 n — 1 


ri 2 n — 1 n — 1 


(A.18) 


and for the derivative of logarithm we obtain: 


da log V = da log 


n 


n -1 


+ ri2qdn 


n -1 
n 


rcoth(r) 11 n^h 
1 + -^77-7 + 


ri 2 n — 1 n — 1 


+ 0(92 
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Finally we can compute the leading order for the V functions: 


= 


2TT‘^q 


1 - ri2q 
r sin(t) 


1 


Va = - 

27r2 2tt‘^ sinh(r) 


r coth(r^ 
ri2 

- T /0 


n -1 

+ 0{q) = V^ + Oiq) 


+ Oiq) 


V3 = 


ri2 n 


+ Oiq) 


27r2 n - 1 

Combining the previous expansion we arrive to the hnal results: 


FVf,Vs = df, ^logn - n2q^ - ri2qdf, + 0{q‘^) 

fv^Va = 2tt\ ((n - i)df,vi + logn) + oiq'^) 

FV^Vs = Oiq^) 


(A.19) 

(A.20) 

(A.21) 

(A.22) 

(A.23) 

(A.24) 


B Details for the analytic form of the Harrington-Shepard determinant 


In this appendix we resume the results obtained in [ 6 ] for the computation of the determi¬ 
nant for caloron gauge fields that satisfy: 

= (B.l) 

We will be interested in the cases a = ru and a = q, with Cq, = ^, ri 2 respectively. 

The variation of the determinant for spin 1 is given by the formula [1, 6 , 24-26]: 


dalogdet {-D"^) = I da' I 4 Aq,/ -|- Ba' + CQ,/-ri2 


where [ 6 ] 


Afy — 

Ba = 
H = 


Ca 1 

12 4^ 

CcK 


dx‘ 


/R3xS1 


ns 


L 


4vr2 Ji?3x5i 

- 1 

h‘^-l + 2f 

h = coth(r)- 

r 

The integration in ri 2 is easily done noting that: 

So we have: 


ns 


n 2 

(1 + ^2/) 

rri 2 

1 

/ 

dr 12 —M2W 

Jo 

^12 

rri 2 

1 

/ 

dr 12—Br^X) 

/o 

^12 


(^ 2-1 + 2 /) 


1 1 ((5Mn)2)' 


12167 r 2 n 4 


1 

167r2 n4 


H{r, t) 


(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B. 6 ) 

(B.7) 

(B. 8 ) 

(B.9) 
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^{ 1 ^ 12 ) contains a logarithmic singularity and can be regularized using the single zero- 
temperature instanton [ 1 , 26] 


A{n2) = 


1 1 
12 16 ^ 


dx 


Ir^xS^ 


n4 


dx 


((5Mno)2)'' 


Ir^ 


TT^ 


(B.IO) 


Numerically, it has been observed [1, 6 ] that the relation B{ri 2 ) = 12 A(ri 2 ) holds, although 
there is not a formal proof of this relation yet. A{ri 2 ) was parametrized in very good 
agreement through the functional forms [Ij: 


^(A) = --3-log fl -h 0.0129 ( 1 -h 0.899 


12 ° V 


6 / 


cfr'^ 


The logarithmic terms dominate over the whole range. 

For the variation in terms of the holonomy (cq = ^ 12 ), we obtain 


dg log det (-D^) |g=o = gri 2 {^A + B + 


where 


A = 

B = 


r 

ns 

ns 


12 47r2 7^3x51 

— f 

4vr2 7^3x51 


(B.ll) 


(B.12) 

(B.13) 

(B.14) 


We find a numerical relation B = 12A — 1. Combining Eqs. (B.9), (4.20), (B.13) and 
(B.14) and using B = ri 2 clri 2 -Sj we obtain our final result for the determinants: 


( 5^12 log det {-Df) = 16 n(ri 2 ) -ri 2 


(B.15) 


6q log det {-D^) |q=o = ^16n(ri2) -h ^ log(ri 2 ) ^ri 2 ^ (B.16) 
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